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THE ENDPOINT REGION IN RADIATIVE QUARKONIA DECAYS ab 
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Department of Physics, Pennsylvania State University, University Park PA 16802 

We consider the inclusive radiative decays quarkonium — > 7 + hadrons and examine 
the effects of soft QCD radiation on the photon energy spectrum near the endpoint. 

The use of perturbation theory to describe radiative decays of quarkonia is based 
on the fact that, as long as the velocity v of the heavy quark is small, these decay 
processes involve two widely separated distance scales: the scale l/(mi> 2 ) over which 
the quark and antiquark bind into the quarkonium and the scale 1/m over which 
the quark-antiquark pair decays (with m the heavy quark mass). By expanding 
about the nonrelativistic limit v — > 0, one may treat the process as the product of a 
long-distance, nonperturbative factor, containing all of the bound-state dynamics, 
and a short-distance factor, describing the annihilation of the heavy quark pair and 
computable as a power series expansion in a s . 

However, near the exclusive boundary of the phase space, where the photon's 
energy E-y approaches its kinematic limit, both the expansion to fixed order in a s 
and the expansion to fixed order in v become inadequate to represent correctly 
the physics of the decay. On one hand, potentially large terms in ln(l — zljwith 
z = Ey/m) appear in the coefficients of the expansion in a s to all orderail. On 
the other hand, classes of relativistic corrections that by power counting are higher 
order in v become enhanced by powers of 1 /a Ju an. These behaviors depend on the 
soft color radiation associated with the decay. In both cases, reliable results can 
only be obtained after resummation of the enhanced contributions. An analogous 
sensitivity to infrared dynamicsis seen in calculations that incorporate models for 
the hadronization of partonsaBQ: nonperturbatiye contributions are found to be 
essential to describe the photon energy spectrumcm near the endpoint. 

In this talk we outline a treatment of the infrared radiation near the kine- 
matic boundary, including soft-gluon coherence effects, and give an argument for 
the cancellation of all the corrections in ln(l — z Xi n the short distance coefficient for 
the color-singlet Fock state in the quarkoniuml!3. We focus on processes in which 
the photon is directly coupled to the heavy quarks, because near the endpoint 
these contributions dominate the contributions in which the photon is produced by 
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fragmentationS of light quarks and gluons. 

Consider the lowest order contribution for a 3 S\ quarkonium H of mass Mm a 
color singlet state, H — > ^gg. The normalized photon spectrum at this order isE3 

J_ d£o _ 1_ rl rl 

Tq dz ir 2 — i 

with 



dx\ I dx2 Mq(xi,X2, z) S(z — 2 + X1 + X2) Q(xx+X2 — 1) , 
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s (1-Xl) 2 (l-X 2 ) 2 {1-zf 

M { X1 ,X2, z) = v 2 + v 2 f + V 2 / . (2) 

Here 2i,a;2 denote gluon energy fractions, Xi — 2Ei/M. The result (|l|) is well 
approximated by a spectrum rising linearly with z. It goes to a constant as z — > 1. 

To higher perturbative orders, logarithmic corrections arise from soft and collinear 
gluon radiation. By power counting the leading behavior of the spectrum as z — > 1 
is of the typea 

1 dT °° 

- const. +y ci a' ln 2fc (l - z) , z -» 1 . (3) 

1 dz f— ' 

fc=i 

The coherent branching algorithmEl provides a method to evaluate the relevant 
multiparton matrix elements in the soft and collinear regions. This algorithm effec- 
tively replaces the calculation of higher-loop Feynman graphs by the calculation of 
tree-level graphs, in which the angular phase space is subject to ordering constraints 
at each branching. The basis for this method is the coherence property of soft gluon 
emission. 

This method has been developed and applied to the study of event-shape ob- 
servables in e + e~ annihilation!!! The photon spectrum in decays of quarkonia can 
be treated in a similar manner. Following this approach, the photon spectrum 
can be expressed, to all orders in a s and to the leading and next-to-leading ac- 
curacy in ln(l — z), in terms of the on-shell amplitude for the tree-level process 
H(P) — > j{k)g(ki)g(k2) andthe mass distributions J g for the time-like jets defined 
by the branching algorithir£3: 



1 dT f d 4 k d A k Y d A k 2 



T dz J (2tt) 3 (2tt) 3 (2tt) 3 



2P-k 



x (2nyS /l (P-k-k 1 -k 2 )S^z- -^j 6+(k z ) 

x M ( - tTee \P,k 1 ,k 2 ) J s ((fci + fc 2 ) 2 ,^) JgHh + kz) 2 ,^) . (4) 

The precise definition of J(p 2 , k 2 ) is given in ref0. The first argument in J is the 
coherence scale; the second argument is the jet mass. To lowest order, J g (p 2 , k 2 ) = 
5{k 2 ) + . . .. In this case in Eq. (0) we reobtain the phase space for the ^ygg final 
state; then the tree-level amplitude M.( tTee > becomes proportional to the amplitude 
A4o of Eq. (^), and Eq. (||) simply gives back the result (0). In general, J satisfies 
an evolution equation of the form 
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J(p 2 ,k 2 )=S(k") + I a a (p l2 )K,(j>'\k 2 )®J(p'\k*) , (5) 

where expressions for the kernel JC are known to leading and next-to-leading orderEl. 
In the present discussion we limit ourselves to considering a double logarithmic 
approximation to J. In this approximationllj 



J dk 2 J g (p 2 , k 2 ) e{Q 2 - k 2 ) « exp 



2 / V 



(6) 



To determine the logarithmic contributions (|^) to the photon spectrum, we need 
to evaluate the branching formula (|]) explicitly. We are interested in the angular- 
ordered, coherent region 

k\ , k 2 < (fci + fc 2 ) 2 < Af 2 . (7) 

The key observation concerns the phase space available for the evolution of the jets 
k\ and k 2 in this region. The boundary on the energy fraction x\ of jet 1 comes 
fromlMa 



i) fragmentation of jet 1: this gives x\ > y/ '4kf/M 2 ; 

ii) recoil of jet 2: after using the angular ordering to approximate the phase 
space in Eq.(^), this gives x\ > k\j\M 2 {\ — z)\. 

For z — * 1 the tightest constraint is set by ii). 

By evaluating the phase space in Eq. (^) explicitly we obtainl!3 



(8) 



1 dT 1 

f ~di ~ 16(27r) 3 (7r 2 -9) 
,i ,i 

dx\ / dx2 Mo(xi,X2, z) 5(z — 2 + xi + x 2 ) <d(xi + x 2 — 1) 
(i Jo 

i 

dfc 2 J g (M 2 (l - z),fc 2 ) e(Af 2 xi(l -z) - k{) &{M 2 xl/A-k{) 

dk 2 J g (M 2 {1 - z), k 2 ) <d(M 2 x 2 (l -z)- 9(M 2 a;^/4 - k\) 

where A4o is given in Eq. (0). 

Eq. (||) allows us to discuss the logarithmic behaviors in the endpoint region by 
using the results for the jet mass distributions. By expanding the double logarithmic 
expression (^) in powers of a s we see that the photon spectrum contains corrections 
involving integrals of the form 



1 7 M 2 (l-z) x 

, . yAftifl-z)/ ,/j ; 



dT i ln ( M % Z \ ) = - I d Xl lnxx . (9) 



That is, logarithmic contributions in x\ arise, which are important at the kinematic 
limit x\ — ► 0, but these never give rise to logarithms of (1 — z) in the spectrum. All 
ln(l — z) cancel because of the form of the constraint on the jet mass, which in turn 
is a consequence of the angular ordering of the gluon radiation. 
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Therefore, higher orders of perturbation theory do not contribute a Sudakov 
suppression of the photon spectrum in the endpoint region. They give rise to a 
constant shift compared to the lowest order answer. 

The picture underlying this result can be understood in simple terms. In the 
boundary kinematics the photon recoils against two almost-collinear gluon jets. The 
cancellation of Sudakov corrections reflects the fact that color is neutralized already 
at the level of this two-jet configuration. This situation may be contrasted with the 
situation one encounters in the decay of an electroweak gauge boson into jetsElfi 
Here Sudakov corrections arise precisely from the presence of color charges in two-jet 
kinematics. 

This picture also indicates that no cancellation should occur for the color-octet 
Fock state in the quarkonium. In this case, we expect the usual Sudakov suppres- 
sion to take place near the endpoint of the photon spectrum. Then one of the 
consequences of the result presented here concerns the ratio of the color-octet to 
the color-singlet contributions. This ratio will be smaller than expected from the 
power counting in v and a s obtained by truncating perturbation theory to fixed 
orderB, owing to the different high order behavior of the perturbation series in the 
two cases. 

To conclude, we observe that resummation formulas of the type can 
be expanded to fixed order in a s and matched with leading and next-to-leadingE-3 
results to obtain improved predictions, valid over a wider range of photon energies. 
If these formulas are combined with models for the infrared behavior of the^irong 
coupling|l3, they can be used to model the nonperturbative shape functionstll that 
parameterize power-like correctionsB'El near the endpoint. Taking account of effects 



from the soft region will likely influcnefflfc thii estimate of the uncertainty on the 
determination of a s from quarkonia decaysaolia. 
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